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As a key to undertanding the basic mechanism for fast reconnection in solar flares, plasmoid- induced- 
reconnection and fractal reconnection are proposed and examined. We first briefly summarize recent solar 
observations that give us hints on the role of plasmoid (flux rope) ejections in flare energy release. We then 
discuss the plasmoid-induced-reconnection model, which is an extention of the classical two-ribbon-flare 
model which we refer to as the CSHKP model. An essential ingredient of the new model is the formation 
and ejection of a plasmoid which play an essential role in the storage of magnetic energy (by inhibiting 
reconnection) and the induction of a strong inflow into reconnection region. Using a simple analytical 
model, we show that the plasmoid ejection and acceleration are closely coupled with the reconnection 
process, leading to a nonlinear instability for the whole dynamics that determines the macroscopic re- 
connection rate uniquely. Next we show that the current sheet tends to have a fractal structure via the 
following process path: tearing =>• sheet thinning => Sweet-Parker sheet =>• secondary tearing =>• further 
sheet thinning => ... These processes occur repeatedly at smaller scales until a microscopic plasma scale 
(either the ion Larmor radius or the ion incrtial length) is reached where anomalous resistivity or colli- 
sionless reconnection can occur. The current sheet eventually has a fractal structure with many plasmoids 
(magnetic islands) of different sizes. When these plasmoids are ejected out of the current sheets, fast 
reconnection occurs at various different scales in a highly time dependent manner. Finally, a scenario is 
presented for fast reconnection in the solar corona on the basis of above plasmoid-induced-reconnection in 
a fractal current sheet. 



1. Introduction 

Recent numerical simulations (e.g., Ugai 1986, 1992, 
Scholer 1989, Biskamp 1986, Yan, Lee, Priest 1992, 
Yokoyama and Shibata 1994, Magara and Shibata 1999, 
Tanuma et al. 1999, 2000) have revealed that if the re- 
sistivity is spatially uniform, fast, steady-state Petscheck- 
type reconnection does not occur but instead slow, Sweet- 
Parker-type reconnection occurs. This holds especially 
when a strong inflow is imposed at the external bound- 
ary, and the only way so far found to achieve a steady- 
state Petschek configuration is to have a localized resis- 
tivity. The so called anomalous resistivity satisfies this 
condition. However, there are a number of questions 
about it. 

1. In order to produce anomalous resistivity, the cur- 
rent sheet thickness must be as small as the ion Larmor 
radius Q 



or the ion inertial length 
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1 Although the physics of anomalous resistivity has not yet 
been fully understood, it is known that anomalous resistivity oc- 
curs due to plasma turbulence which is produced by the micro- 
scopic plasma instability, such as the lower hybrid drift instability, 
the electrostatic ion cyclotron instability, and the ion sound in- 
stability (e.g., Treumann and Baumjohann 1997). In the case of 
the lower hybrid drift instability, the threshold of the instability 
is Vd > Ujonth where v^ = j/(n e e) is the electron-ion drift ve- 
locity, and v ion t h = (fcT/m;) 1 / 2 is the ion thermal speed. This 
equation becomes equivalent to d < Tjj f on if we consider the 
pressure balance p = InkT ~ B 2 /Sir between inside and outside 
of the current sheet, where d is the thickness of the current sheet. 
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both of which are of order of 1 m in the solar corona. 
Since the size of solar flares is typically 10 4 km, there is 
a large gap between the flare size and the necessary mi- 
croscopic scale to produce anomalous resistivity. How 
can such an enormous gap between macroscopic and 
microscopic scales be reconciled in real flares? 

2. Even if the anomalous resistivity (or localized re- 
sistivity) is realized, what determines the reconnection 
rate? 

Based on recent observations of solar flares and nu- 
merical simulations, we try to give possible answers to 
above questions. We argue that the key physics needed 
to answer the above questions is the global coupling 
between plasmoid (flux rope) ejection and reconnection 
process. Since this coupling is scale free, it can occur 
on any scale, constituting a fractal reconnection pro- 
cess, which couples the macro- and micro-scales. 
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2. Solar Observations: Flares and Plas- 
moid Ejections 

Yohkoh has revealed numerous indications of mag- 
netic reconnection in solar flares, such as cusps, arcades, 
loop top hard X-ray (HXR) sources, X-ray jets and so 
on (e.g., Tsuneta et al. 1992a, Hanaoka et al. 1994, 
Masuda et al. 1994, Forbes and Acton 1996, Shibata 
1999). Furthermore, as has been predicted by some pi- 
oneers (Hirayama 1991, Moore and Roumcliotis 1992), 
the association of plasmoid (flux rope) ejections with 
flares is much more common than had been thought 
(e.g., Shibata et al. 1995, Nitta 1996, Ohyama and 
Shibata 1997, 1998, 2000, Tsuneta 1997, Akiyama and 
Hara 2000). This has led us to advocate a unified 
model of flares shown in Figure 1 (Shibata et al. 1995, 
Shibata 1996, 1997, 1998, 1999). Recent observations 
with SOHO/LASCO have also revealed a lot of evi- 
dence of flux rope and disconnection events in coronal 
mass ejections (CMEs) (e.g., Dere et al. 1999, Simnet 
et al. 1997), and Yohkoh has shown that giant arcades 
formed after prominence eruptions or CMEs are phys- 
ically similar to flare arcades even though their total 
X-ray intensity is much lower than that of normal flares 
(e.g., Tsuneta et al. 1992b, Hanaoka et al. 1994). Fig- 
ure 2 shows several examples of plasmoid (flux rope) 
ejections on the Sun from the largest scale in CMEs 
(<~ 10 11 cm) to the smallest scale in compact flares 
(<~ 10 9 cm). The velocity of these plasmoids range from 
a few 10 km/s to 1000 km/s, and their maximum val- 
ues are comparable to the inferred coronal Alfven speed 
(~ 1000 km/s). These images show that the magnetic 
reconnection and associated plasmoid ejection univer- 
sally occur on widely different scales. 

One of the interesting findings by Yohkoh concerning 
X-ray plasmoid ejections is that, in impulsive flares, 
a plasmoid starts to be ejected slowly, long before the 
impulsive phase, and then is rapidly accelerated during 
the impulsive phase (Ohyama and Shibata 1997, 1998, 
Tsuneta 1997; Fig. 3). Similar behavior has also been 
observed for LDE flares and CME events (e.g., Kahler 
et al. 1988, Hundhausen 1999). 

Another interesting finding from Yohkoh on X-ray 
plasmoid ejection is that there is a positive correlation 
between the plasmoid velocity (V p i a smoid ~ 30 — 400 
km/s) and the apparent rise velocity of the flare loop 
(Vioop ~ 4 - 20 km/s) (Shibata et al. 1995): 

Vplasmoid — 

(8 - 20) x V loop . (3) 

This relation (though still very preliminary) suggests 
that the plasmoid velocity is related to the reconnec- 
tion inflow speed, or vice-versa. This is because the 
apparent rise motion of the flare loop is coupled to the 
reconnection process. Consequently, magnetic flux con- 
servation leads to 

Vloop — {-^inflow I Bloop)Vi n f low • (4) 

Morimoto and Kurokawa (2000) found a correlation be- 
tween the erupting velocity of H-alpha filaments (i.e., 
a plasmoid) and the thermal energy density of post- 
eruption X-ray arcades. This also suggests that there 



is a physical relation between plasmoid velocity and in- 
flow speed (reconnection rate). 

3. Role of Plasmoid: Plasmoid-Induced- 
Reconnection Model 

On the basis of these observations, Shibata (1996, 
1997) proposed a plasmoid-induced-reconnection model, 
which is an extension of the classical CSHKP (Carmichcl 
1964, Sturrock 1966, Hirayama 1974, Kopp and Pneu- 
man 1976) model and similar in spirit to the model of 
Anzer and Pneuman (1982). In this model, the plas- 
moid ejection plays a key role in triggering fast recon- 
nection in two different ways (Fig. 1). 

1) A plasmoid (flux rope) can store energy by inhibit- 
ing reconnection. A large magnetic island (plasmoid or 
flux rope) inside the current sheet is a big obstacle for 
reconnection. Hence if an external force compresses the 
current sheet, magnetic energy can be stored around 
the current sheet. Only after the plasmoid is ejected 
out of the current sheet, will the anti-parallel field lines 
be able to touch and reconnect. If a larger plasmoid is 
ejected, a larger energy release occurs. 

2) A plasmoid ejection can induce a strong inflow 
into the reconnection site. If a plasmoid is suddenly 
ejected out of the current sheet at the velocity V p i aS moid, 
an inflow must develop toward the X-point in order to 
compensate for the mass ejected by the plasmoid, as has 
been shown in many numerical simulations (e.g., Forbes 
1990, Yokoyama and Shibata 1994, 2000, Magara et al. 
1997, Tanuma et al. 2000; see also Fig. 4). The inflow 
speed can be estimated from the mass conservation law 
(assuming incompressibility, for simplicity); 

Vinflow ^ Vpl asrno idWplasmoid I ^inflow, (5) 

where W p i aS moid is the typical width of the plasmoid, 
and L in fi ow (> W p i asmoi d) is the typical vertical length 
of the inflow region. In deriving equation (5), it is 
assumed that the mass flux into reconnection region 
(~ Li n fi ow Vinfi ow ) is balanced by the mass flux carried 
by the plasmoid motion (~ V p i asmo idW p i asmold ). Since 
the reconnection rate is determined by the inflow speed, 
the ultimate origin of fast reconnection in this model is 
the fast ejection of the plasmoid. If the plasmoid ejec- 
tion (or outflow) is inhibited in some way, then fast 
reconnection ceases (Ugai 1982, Tanuma et al. 2000, 
Lin and Forbes 2000). 

This model naturally explains (1) the strong accel- 
eration of plasmoids during the impulsive (rise) phase 
of flares (see Fig. 3 and next section), (2) the positive 
correlation between plasmoid velocity and the apparent 
rise velocity of flare loops (eqs. 3 and 5), (3) the total 
energy release rate of flares and plasmoid ejections (Shi- 
bata 1997), and (4) the time scale of the impulsive (rise) 
phase for both impulsive flares (~ L in fi ow /V p i asmold ~ 
10 4 km/100 km/s - 100 sec), and for LDE flares (~ 10 5 
km/100 km/s - 10 3 sec). 

It is interesting to note that similar impulsive recon- 
nection associated with plasmoid ejection (current sheet 
ejection) has also been observed in laboratory experi- 
ments (Ono et al. 2000). 
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4. Nonlinear Instability caused by Strong 
Coupling between Plasmoid Ejection 
(Acceleration) and Reconnection 

In this section, we examine the physical mechanism of 
the plasmoid-induced-reconnection in more detail. We 
consider a situation where reconnection has just begun 
and a plasmoid, with a length L p and a width W p , has 
just started to form. The reconnection generates a jet 
(with the Alfven speed Va) which collides with the plas- 
moid and accelerates it. Thus the plasmoid speed in- 
creases with time, which induces a faster inflow into the 
reconnection point (i.e. the X-point), thereby leading 
to yet faster reconnection and an even larger energy re- 
lease rate. This, in turn, accelerates the plasmoid again, 
eventually leading to a kind of nonlinear instability for 
the plasmoid ejection and the associated reconnection. 

Let us estimate the plasmoid velocity in this process, 
by assuming that the plasmoid is accelerated solely by 
the momentum of the reconnection jet. (Note that we 
do not deny the possibility of acceleration of plasmoid 
by other mechanism such as global magnetic pressure. 
The purpose of this section is simply to show how the 
momentum of the reconnection jet can accelerate the 
plasmoid.) We also assume that the plasmoid den- 
sity p p and the ambient plasma density p are constant 
with time, for simplicity. In absence of any appropriate 
time-dependent theory in a rapidly evolving configura- 
tion, we assume that the steady state mass conservation 
ViLi = VpWp (equation (5)) is valid and also that all 
the mass flux (ViLi) convected into the reconnection 
region (with length of Li) are accelerated up to Alfven 
speed as in Sweet-Parker or Petschek model. 

We first consider the case in which the mass added 
to the plasmoid by the reconnection jet is much smaller 
than the total mass of the plasmoid (i.e., the plasmoid 
speed V p is much smaller than the Alfven speed Va)- 
Equating the momentum added by the reconnection jet 
with the change of momentum of the plasmoid, we have 

PpLpWp^- = pVL 2 V A = pVpWpV A (6) 

where we use the mass conservation relation for the 
inflow and the outflow, V P W P = ViLi (eq. 5). (See 
Appendix for detailed derivation of the equation (6).) 
Physically, this means that the inflow is induced by the 
outflow (plasmoid ejection). This is the reason why this 
reconnection is called plasmoid-induced-reconnection. 
The equation (6) is easily solved to yield the solution 

V p = Va exp(^) (7) 

where Vo is the initial velocity of the plasmoid, and 

pV A 



PpLp 



(8) 



Thus, the plasmoid velocity increases exponentially with 
time, and the "growth time" (1/lu) is basically of order 
of Alfven time. The inflow speed becomes 



y = Wp = W p V exp(ut) 
1 Li p L l (0) + ^(expM)-l) 



If Wp is constant, the inflow speed increases exponen- 
tially with time in the initial phase, but tends to be a 
constant (~ ujW p ) in the late phase. 

As time goes on, the mass added to the plasmoid by 
the jet increases and eventually becomes non-negligible 
compared with the initial mass (i.e., the plasmoid speed 
becomes non- negligible compared with the Alfven speed) . 
In this case, we obtain the solution (see Appendix for 
derivation) : 



_ Va exp(o;f ) 
p exp(wt) — 1 



Va/Vq ' 



(10) 



Hence the plasmoid speed is saturated at around t = 
t c ~ ^\ii(Va/Vo) and hereafter tends to the Alfven 
speed Va as time goes on. The inflow speed becomes 

y WpVp 



Va exp(wt) / (exp(u>t) + a) 



(11) 



p {Va/u) ln[(exp(wt) + a)/(l + a)] + L,(0) 

where a = Va/Vq — 1. If W p is constant in time, the 
inflow speed gradually decreases in proportion to 1/i 
after t c . ^ On the other hand, if W p increases with time 
in proportion to t after t Cl the inflow speed becomes 
constant, 



V = L0Wp{t = 0) 



pV A 
PpLp 



W p (t = 0) 



(12) 



In this case, the reconnection becomes steady, and the 
shape of the reconnection jet and plasmoid becomes 
self-similar in time and space (e.g., Nitta et al. 2000, 
Yokoyama and Shibata 2000). 

A typical solution for W p = constant is shown in 
Figure 5, which reminds us of the observed relation be- 
tween plasmoid height vs. hard X-ray intensity (Fig. 3; 
Ohyama and Shibata 1997) and explains also the nu- 
merical simulation results (Fig. 4; Magara et al. 1997) 
very well. It is noted here that the hard X-ray intensity 
is a measure of either the electric field at the reconnec- 
tion point (E oc V.B) or the energy release rate (oc 
Poynting flux oc VB 2 /(4tt)). 

5. Fractal Reconnection 

As we discussed in section 1, we have a fundamental 
question: how can we reach the small dissipation scale 
necessary for anomalous resistivity or collisionless re- 
connection in solar flares? Also, even if we can reach 
such a small scale, is it true that there is only one dif- 
fusion region with a thickness of 100 cm (and with a 
length of 10 m or 100 m) in a solar flare as expected 
from Petschek's steady state theory? 

The idea that the reconnection process is inherently 
turbulent, involving a spectrum of different scales, has 
been around for some time (see Ichimaru 1975, for ex- 
amples). However, here we argue that a plasma with 



(9) 



This kind of evolution occurs when 1) the current sheet length 
is limited (Tanuma et al. 2000), 2) magnetic field distribution is 
non-uniform around the current sheet (Magara et al. 1997). 
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large magnetic Reynolds number (occurs as in the so- 
lar corona, the interstellar medium, or the intergalac- 
tic medium) inevitably leads to a fractal current sheet 
with many magnetic islands of different sizes connecting 
macroscopic and microscopic scales (Tajima and Shi- 
bata 1997, Shibata et al. 1997, 1998, Tanuma et al. 
2000). 

Let us first consider the Sweet-Parker current sheet 
with a thickness of S n and a length A„. This current 
sheet becomes unstable to secondary tearing if 



t n < A„/Va, 



(13) 



where t n is the growth time of the tearing instability at 
maximum rate (iOmax ^ kmax 2/5 tdif~ 3/5 t A ~ 2/5 and 
kmax oc {t dlf /t A )~ 1/4 , where ^max and k max are the 
maximum growth rate and corresponding wave num- 
ber), 

/f) 2 /) \ 1/2 

t n x (t dlf t A )^ * , (14) 

and Xu/Va is the time for the reconnection flow to carry 
the perturbation out of the current sheet. (As for the 
theory of the secondary tearing in the Sweet-Parker 
sheet, see e.g., Sonnerup and Sakai (1981), Biskamp 
(1992).) That is, if t n > Xu/Va, the tearing instability 
is stabilized by the effect of flow. Using eq (13) and 
(14), we find 

'A„\ 2 



i.e., 



(15) 



If this inequality is satisfied, the secondary tearing oc- 
curs, leading to the current sheet thinning in the non- 
linear stage of the tearing instability. At this stage, the 
current sheet thickness is determined by the most un- 
stable wavelength of the secondary tearing instability, 
i.e., 

X n+1 ~ 66 n R£, n = 6^/4^/^5/4 < 6r ? 1 /6y-V6 A 5/6 5 

(16) 

where R m *,n — 8 n V A /r). The current sheet becomes 
thinner and thinner, and when the current sheet thick- 
ness becomes 



(17) 



further secondary tearing occurs, and the same process 
occurs again at a smaller scale (Fig. 6). It follows from 
eq (16) and (17) that 



or 



where 



and 



*^(i)">«. 

s -i <- <¥r 

A = 6 2 / 3 R m 1 / 6 , 
LV A 



R„ 



V 



(18) 

(19) 
(20) 
(21) 



This fractal process continues until the current sheet 
thickness reaches the microscopic scale such as the ion 
Larmor radius or ion inertial length. The equation (19) 
leads to 

^ = A 6(l-*)(|)* 5 (22a) 



where 



x = (5/6)". 



(226) 



From this, we can estimate how many secondary tear- 
ings are necessary for the initial macroscopic current 
sheet to reach the microscopic scale. Taking the typ- 
ical solar coronal values, So = 10 8 cm, L = 10 9 cm, 
V A = 10 8 cm/s, r] = 10 4 cm 2 /s for T = 10 6 K, wc find 
R m = 10 13 and 

A ~ 0.02. (23) 

Since S n must be smaller than the typical microscopic 
scale, e.g., the ion Larmor radius (~ 100 cm), we have 



S n /L < 



(24) 



or 



(0.02) 6 < 1 -( 5 / 6 )")(0.1)( 5 / 6 )" < 10- 7 . 
The solution of this inequality (see Fig. 7) is 

n > 6. (25) 

That is, in the solar corona, six secondary tearings are 
necessary to reach microscopic current sheet. 

What is the time scale of this fractal tearing? The 
time scale for the n-th tearing is 

t n ~ ST^Va)- 1 ' 2 = (S n /S ) 3/2 to, (26) 



(27) 



where 

to = 6 3 / 2 /(r,V A ) 1 / 2 . 
Since equation (22) leads to 

5 n /6 *4 (1 - (5/6r \ (28) 
where A = G 2 ^ 3 R m 1 J_Q, and i? m » ;0 = SoVa/t], we find 

t n ~ ^ (1 - (5/6) " ) t . (29) 

Thus we obtain 

W*n-1 ^ 4 3/2)(5/6) "" < AT (30) 

for n > 1. It follows from this equation that 



tn ^ A J t n -i < A, 



(3/2), 



(31) 



Consequently, the total time from the 1st (secondary) 
tearing (t\) to the n-th (secondary) tearing (t n ) be- 
comes 



ttotai =h+t 2 + ... + t n < t A 3 a /2 1 A ° < t A 3 /2 . 



3n/2 



(32) 

For typical coronal conditions (described above), this 
time scale becomes 



ttotai < 6 x 10-%, 



(33) 
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which is much shorter than the time scale of the 0-th 
tearing (to)- Although the 0-th tearing time is long (~ 
3 x 10 4 — 10 6 sec for initial current sheet with So ~ 10 7 — 
10 8 cm), the nonlinear fractal tearing time is quite short 
(less than 10 2 — 3 x 10 3 sec), so that the microscopic 
scale is easily reached within a short time as a result of 
the fractal tearing. 

It should be stressed that the role of the fractal tear- 
ing is only to produce a very thin current sheet with 
a microscopic scale of order of the ion Larmor radius 
or the ion inertial length. The fractal tearing does not 
explain the main energy release in flares. The main en- 
ergy release is explained by the fast reconnection pro- 
cess which occurs after the ejection of the large scale 
plasmoid as we discussed before. 

6. Summary : A Scenario for Fast Re- 
connection 

Let us summarize our scenario of fast reconnection 
in the solar corona, which is illustrated in Fig. 8 and 
Fig. 9 (the latter is from a numerical simulation by 
Tanuma et al. 2000 and it nicely illustrates a part of our 
scenario) . Our scenario can also be applied to other hot 
astrophysical plasmas (e.g., stellar corona, interstellar 
medium, galactic halo, galactic clusters, and so on) for 
which magnetic Reynolds number and the ratio of its 
characteristic scale length to the ion Larmor radius (or 
ion inertial length) are very large. 

Initially we assume the current sheet whose thick- 
ness is much larger than the microscopic plasma scale. 
Such a current sheet is easily created by the interaction 
of emerging flux with an overlying coronal field (e.g., 
Heyvaerts et al. 1977, Shibata et al. 1992, Yokoyama 
and Shibata 1995), the collision of a moving bipole with 
other magnetic structure (e.g., Priest et al. 1994), the 
global resistive MHD instability in a shearing arcade 
(e.g., Mikic et al. 1988, Biskamp and Welter 1989, Ku- 
sano et al. 1995, Choe and Lee 1996, Magara et al. 
1997, Choe and Chen 2000), or other related mech- 
anisms (e.g., Forbes 1990, Chen, Shibata, Yokoyama 
2000). 

If the current sheet length becomes longer than the 
critical wavelength for the tearing mode instability, the 
instability starts. As the instability developes, it enters 
a nonlinear regime which makes the initial current sheet 
thinner and thinner. The current sheet thinning stops 
when the sheet thickness becomes comparable to that of 
the Sweet-Parker sheet, and thereafter the sheet length 
increases with time. If the sheet length becomes longer 
than a critical wavelength (eq. 13), secondary tearing 
occurs. Even if the sheet has not yet reached the Sweet- 
Parker state, it can become unstable to the secondary 
tearing if the sheet thickness satisfies the same condi- 
tion (eq. 13). Then the same process occurs again at 
a smaller scale, and the system evolves into one that is 
fractally structured. In this way, a microscopicly small 
scale (such as ion Larmor radius or ion inertial length) 
can be reached within a short time. 

Once a small scale is achieved, fast reconnection oc- 
curs because anomalous resistivity can now set in. It is 



also possible that fast collisionless reconnection occurs 
with a nondimcnsional reconnection rate of the order 
of 0.1-0.01 at this small scale (see recent full particle 
simulations by, e.g., Drake 2000, Hoshino et al. 2000, 
Tanaka 2000, Horiuchi and Sato 2000). Hence small 
scale magnetic islands (plasmoids) created by small scale 
tearing are ejected at the Alfven speed and collide with 
other islands to coalesce with each other, thereby mak- 
ing bigger islands (plasmoids). This coalescing process 
itself also occurs with a fractal nature (Tajima and Shi- 
bata 1997). 

It should be noted that the ejection (acceleration) of 
plasmoids (flux rope with axial field in 3D space) can 
enhance the inflow into the reconnection point, creat- 
ing a positive feedback, i.e., nonlinear instability (as we 
outlined in section 4). This determines the macroscopic 
reconnection rate which may be smaller or larger than 
the microscopic reconnection rate. If the macroscopic 
reconnection rate (inflow speed) is larger than the mi- 
croscopic reconnection rate, the magnetic flux is accu- 
mulated around the diffusion region, leading to inter- 
mittent fast reconnection (Lee and Fu 1986, Kitabata 
et al. 1996, Schumacher and Klicm 1996, Tanuma et al. 
1999, 2000; Fig. 9). On ther other hand, if the macro- 
scopic reconnection rate is smaller than the microscopic 
reconnection rate, the reconnection may continue in a 
quasi-steady state. However, there may be large am- 
plitude perturbations around the reconnection point, 
so that it would be difficult to maintain quasi-steady 
reconnection. The reconnection would be very time de- 
pendent with intermittent reconnection and ejection of 
plasmoids with various sizes created by fractal recon- 
nection. Petschek's slow shocks are also formed in a 
very time dependent manner (e.g., Yokoyama and Shi- 
bata 1994, Tanuma et al. 2000). The local macroscopic 
reconnection rate can be much larger than the average 
reconnection rate and is determined by the macroscale 
dynamics, i.e., plasmoid-induced-reconnection. In this 
case, the time dependence is essential for determining 
the reconnection rate. 

Since this process is scale free, we have fractal struc- 
ture in the global current sheet. The greatest energy re- 
lease occurs when the largest plasmoid is ejected. This 
may correspond to the impulsive phase of flares. The 
time variation of the reconnection rate (and the total 
energy release rate) associated with ejection of plas- 
moids with various sizes is also fractal. That is, the 
power spectrum of the time variation of the reconnec- 
tion rate and the energy release rate show a power-law 
distribution. This may correspond to the fragmented 
light curves of solar X-ray and radio emissions in the 
impulsive phase of flares (e.g., Benz and Aschwandcn 
1992). 

Quantitative proof of the fractal nature of the cur- 
rent sheet (especially in 3D geometry) remains as an 
important subject for future numerical simulations and 
laboratory experiments, both of which will have to be 
able to handle much larger magnetic Reynolds number 
than they currently do (i.e. R m ~ 10 3 — 10 4 ) in order 
to solve this fundamental problem. 
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Appendix: Derivation of Equations (6) 
and (10) 

As we wrote in the text, we assume that all the mass 
convected into the reconnection region (pViLi per unit 
time per unit length in 2D space) are accelerated to 
Alfven speed Va- Since such accelerated mass (recon- 
nection jet) collides with the plasmoid, it can accelerate 
the plasmoid. Denoting V v — plasmoid speed, M p = 
plasmoid mass, AM p = mass convected by the recon- 
nection jet during a short time At, which is equal to 
increase in plasmoid mass during At, we obtain the 
conservation of momentum as 

AM P V A + M P V P = (M p + AM P )(V P + AV P ). (Al) 

Here the left hand side is the total momentum before 
collision, and the right hand side is the total momentum 
after collision. If we neglect the term AM p in the right 
hand side of equation (Al) (i.e., if we assume V p <C Va), 
we have 

M P AV P = AM P V A . (A2) 

The plasmoid mass (M p ) and the mass added to the 
plasmoid (AM p ) by the jet for a short time At are 
written as 

M p = PpLpWp, (A3) 

AMp = pViLiAt, (AA) 

both of which are per unit length. Using these formulae, 
the equation (A2) becomes p p L p WpAVp = pViLiVAAt, 
which is equivalent to 

PpLpWp^ = P V l L l V A . (Ah) 

Since we assumed that the mass injection into re- 
connection region is induced by the plasmoid motion, 
i.e., ViLi = VpWp (equation (5)), the right hand side of 
equation (A5) becomes equal to pV p W p Va, so that we 
get equation (6): 

PpLpWp^- = pViLiV A = pVpWpV A . (6) 

Note that in deriving above formulae, we did not as- 
sume conservation of kinetic energy. This is because 
some part of the kinetic energy is dissipated to heat 



the plasmoid, leading to increase in gas pressure (inter- 
nal energy) of the plasmoid. Although such enhanced 
gas pressure may accelerate the plasmoid further, we 
neglected the effect of gas pressure in above treatment 
for simplicity, since it is not easy to estimate how much 
fraction of internal energy is converted to the kinetic 
energy of a plasmoid. 

When V p grows to the value that cannot be neglected 
compared with Va, we cannot neglect the term AM p V p 
in the right hand side of equation (Al). In this case, the 
momentum conservation equation becomes M p AV p = 
(Va — Vp)AM p . Combining this equation with equations 
(A3) and (A4), we get 

If p,p p ,Va, and L p are constant in time, the solution 
of this equation becomes equation (10): 

V = VAeMut) (1Q1 

p cxp(ujt)-l + V A /Vo' ^ ' 

Here Vq is the initial velocity of a plasmoid at t = 0, 
and u) = (p/pp)(V A /Lp). 
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Figure Captions 

Fig. 1 A unified model of flares: plasmoid-induced- 
reconnection model (Shibata ct al. 1995). This is an 
extention of a classical model of eruptive solar flares, 
called the CSHKP model. 

Fig. 2 Various plasmoids (flux rope) with different 
scales observed on the Sun. (a) Coronal mass ejection 
(CME), the largest-scale plasmoid on the Sun (~ fO 11 
cm) observed with SOHO/LASCO on Nov. 1-2, 1997 
(Dere et al. f999). These are running-difference im- 
ages. The velocity of the CME is 140 - 240 km/s. 
(b) Large-scale X-ray plasmoid associated with an LDE 
(long duration event) flare on Feb. 21, 1992 (~ 10 10 cm) 
observed with Yohkoh/SXT (Hudson 1994, Ohyama 
and Shibata 1998). The plasmoid velocity is about 100 
km/s. (c) Small-scale X-ray plasmoid associated with 
an impulsive flare (~ 10 9 cm) observed with Yohkoh/SXT 
on Oct. 5, 1992 (Ohyama and Shibata 1998). These are 
negative images. The velocity of the plasmoid is 250 - 
500 km/s. 1" corresponds to 726 km. 

Fig. 3 Temporal variations of the height of an X-ray 
plasmoid and the hard X-ray intensity of an impulsive 
solar flare on 11 Nov. 1993 observed with Yohkoh SXT 
and HXT (Ohyama and Shibata 1997). 

Fig. 4 Temporal variations of both the reconnec- 
tion rate (electric field at the reconnection-point) and 
the height of the plasmoid (magnetic island) for a typ- 
ical result of 2.5D MHD numerical simulation of mag- 
netic reconnection induced by plasmoid (flux rope) ejec- 
tion (Magara, Shibata, Yokoyama 1997). Units of the 
height, time, and electric field are L (a half length be- 
tween footpoints of a sheared arcade loop), to = L/C s q 
(C s o is the sound speed ~ 0.4Va), and E = C s qB /c, 
respectively. In a typical solar coronal condition, L ~ 
5000 km, V A ^ 1000 km/s, E a ~2x 10 4 V/m, and t ~ 
20 sec. 

Fig. 5 Temporal variations of the plasmoid velocity 
(Vpiasmmd), its height, and inflow velocity (V ln fi ow ), in 
an analytical model (eqs. 10 and 11) for the case of 
Va/Vo — 100. Units of the velocity, height, and time 
are Va,L p and L p /Va, respectively. 

Fig. 6 Schematic view of fractal reconnection. 

Fig. 7 The current sheet thickness (S n /L) in the n-th 
secondary tearing (see eq. 22a). 

Fig. 8 A scenario for fast reconnection. I: The ini- 
tial current sheet. II: The current sheet thinning in the 
nonlinear stage of the tearing instability or global resis- 
tive MHD instability. The current sheet thinning stops 
when the sheet evolves to the Sweet-Parker sheet. Ill: 
The secondary tearing in the Sweet-Parker sheet. The 
current sheet becomes fractal because of further sec- 
ondary tearing as shown in Fig. 6. IV: The magnetic 



islands coalesce with each other to form bigger mag- 
netic islands. The coalescence itself proceeds with a 
fractal nature. In the phases III and IV, the micro- 
scopic plasma scale (ion Larmor radius or ion inertial 
length) is reached, so that fast reconnection becomes 
possible at small scales, V: The greatest energy release 
occurs when the largest plasmoid (magnetic island or 
flux rope) is ejected. The maximum inflow speed (Vj n 
= reconnection rate) is determined by the velocity of 
the plasmoid (V p ). Hence this reconnection is called 
plasmoid-induced-reconnection. 

Fig. 9 Numerical simulation of reconncction trig- 
gered by an MHD fast mode shock (Tanuma et al. 
2000), illustrating a part of our proposed scenario for 
fast reconnection (Fig. 8): a) passage of the MHD fast 
shock, b) current sheet thinning (in the nonlinear stage 
of the tearing instability), c) the Sweet-Parker recon- 
nection, d) secondary tearing, e) Petschck reconncction 
as a result of the onset of anomalous resistivity. Slow 
shocks inherent to Petschek reconnection are formed. 
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